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47 Abstract 

48 The coplane analysis technique for mapping the three-dimensional wind field of 

49 precipitating systems is applied to the NASA High Altitude Wind and Rain Airborne Profiler 

50 (HIWRAP). HIWRAP is a dual-frequency Doppler radar system with two downward pointing 

51 and conically scanning beams. The coplane technique interpolates radar measurements to a 

52 natural coordinate frame, directly solves for two wind components, and integrates the mass 

53 continuity equation to retrieve the unobserved third wind component. This technique is tested 

54 using a model simulation of a hurricane and compared to a global optimization retrieval. The 

55 coplane method produced lower errors for the cross-track and vertical wind components, while 

56 the global optimization method produced lower errors for the along-track wind component. 

57 Cross-track and vertical wind errors were dependent upon the accuracy of the estimated 

58 boundary condition winds near the surface and at nadir, which were derived by making certain 

59 assumptions about the vertical velocity field. The coplane technique was then applied 

60 successfully to HIWRAP observations of Hurricane Ingrid (2013). Unlike the global 

61 optimization method, the coplane analysis allows for a transparent connection between the radar 

62 observations and specific analysis results. With this ability, small-scale features can be analyzed 

63 more adequately and erroneous radar measurements can be identified more easily. 
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70 1 . Introduction 

7 1 The use of airborne Doppler radars has significantly advanced our understanding of 

72 meteorological phenomena by providing wind structure information that details the dynamics of 

73 an evolving system. Airborne platforms have been particularly important for observing 

74 phenomena that occur in remote areas, such as tropical cyclones over the open ocean. Lhennitte 

75 (1971) first discussed the idea of using airborne Doppler radars for obtaining three-dimensional 

76 wind structures. A single Doppler radar beam measures the along-beam velocity component of 

77 precipitation particles within that beam. In order to retrieve all three components of the wind 

78 field, two (or more) Doppler radar beams must scan an area with a sufficient separation angle 

79 between the beams (e.g. Armijo 1969; Klimowski and Marwitz 1992). Airborne radars must 

80 therefore employ certain scanning techniques that provide multiple views of the wind from 

81 sufficiently different angles in order to map the wind structure of the precipitation phenomena. 

82 In one of the first airborne Doppler studies, Marks and Houze (1984) utilized a scanning 

83 technique for successful mapping of the three-dimensional wind field. They used data collected 

84 by the X-band Doppler radar on board the National Oceanic and Atmospheric Administration’s 

85 (NOAA) WP-3D (P3) aircraft. Located in the tail of the aircraft, the radar antenna pointed 

86 orthogonally to the aircraft track and scanned circularly through all elevation angles around a 

87 horizontal axis (for more on this radar, see Jorgensen 1984). Multiple viewing angles of the same 

88 domain were obtained by flying the aircraft at different track angles. Another tail radar was later 

89 installed on the second NOAA P3 aircraft allowing for simultaneous Doppler observations when 

90 both aircraft were flown together (Gamache at al. 1995). Both tail radars soon implemented the 

91 fore/aft scanning technique (FAST; Jorgensen and DuGranrut 1991), in which the antenna 

92 alternately points -20° to the fore and aft of the aircraft while circularly sweeping around a 
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93 horizontal axis. With this technique, multiple along-beam velocity measurements from the same 

94 domain are obtained along a single flight track by the different fore and aft angles. The National 

95 Center for Atmospheric Research (NCAR) Electra Doppler Radar (ELDORA) operates with the 

96 same scanning geometry but utilizes two antennas that rotate at a faster rate, allowing for higher 

97 resolution observations (Hildebrand et al. 1996). 

98 Multiple techniques for retrieving the three-dimensional wind field have been developed 

99 for the NOAA P3 tail radar and the ELDORA radar, which both scan around a horizontal axis. 

100 The first method is a local solver that interpolates radial velocities from each viewing angle to a 

101 Cartesian grid and solves for the corresponding velocities in the horizontal plane. These 

102 horizontal velocities from different viewing angles are then used to calculate two orthogonal 

103 horizontal wind components (e.g. Jorgensen et al. 1983; Marks and Houze 1984). The vertical 

104 wind component is calculated by integrating the anelastic mass continuity equation using 

105 appropriate boundary conditions. This technique is simple and computationally inexpensive, but 

106 errors can accumulate in the wind component along the direction of integration (Gao et al. 1999). 

107 The second method is a global optimization approach that minimizes a cost function 

108 containing the differences between the radar-measured and retrieved velocity components. This 

109 cost function also includes constraints such as the anelastic mass continuity equation and vertical 

110 velocity boundary conditions (Gamache 1997; Bousquet and Chong 1998; Reasor et al. 2009). 

111 With the avoidance of explicit integration, this variational technique reduces errors in the vertical 

112 velocity for the aforementioned scanning geometry (Gao et al. 1999). Since all retrieval 

113 strategies are limited by the geometry of the scanning technique, no individual retrieval method 

114 is perfect; however, utilizing multiple methods adds to the reliability of scientific interpretations 

115 of retrieved wind fields. 
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116 The High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP), recently 

117 developed at NASA Goddard Space Flight Center, is a Doppler radar system that employs a 

118 different scanning strategy from the previously mentioned airborne Doppler radars (Li et al. 

119 2014). It operates with two beams that point downward at fixed angles (30° and 40° away from 

120 nadir) with each beam scanning conically around a vertical axis. HIWRAP flew for the first time 

121 in 2010 on the NASA Global Hawk unmanned aircraft during the Genesis and Rapid 

122 Intensification Processes (GRIP) field experiment (Braun et al. 2013). 

123 Recent studies have begun exploring how established retrieval methods can be applied to 

124 the scanning geometry of HIWRAP. Tian et al. (2015) applied the Velocity Azimuth Display 

125 (VAD) technique (Lhermitte and Atlas 1961; Browning and Wexler 1968) to HIWRAP data to 

126 obtain the mean vertical profile of the horizontal wind along the flight track. Under the 

127 assumptions that the wind field is linear and the vertical velocity is constant across the scan 

128 circle, this method fits the measured radial winds at each altitude to a sinusoidal curve as a 

129 function of azimuth. Guimond et al. (2014) implemented the global optimization technique to 

130 obtain the three-dimensional wind field in the HIWRAP scanning domain. The cost function for 

131 this variational scheme included a modified weighting parameter that was better suited for the 

132 different scanning geometry. 

133 In this paper, we extend the application of established retrieval techniques to the 

134 HIWRAP geometry by focusing on a simple interpolation and integration approach. The vertical 

135 integration scheme used for the P3 tail radars cannot be applied in the HIWRAP case since the 

136 scanning geometry does not align sufficiently with the horizontal wind in order to avoid large 

137 projection errors introduced by the vertical wind. A better alternative is the coplane method 

138 described by Armijo (1969) and Miller and Strauch (1974). The coplane method uses a 
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139 cylindrical coordinate system in which two components of the wind are readily derived from the 

140 observations. The third wind component is completely unobserved by the radar and must be 

141 retrieved by explicitly integrating the mass continuity equation with specified boundary 

142 conditions. As a local solver, the solution of the coplane method at a certain gridpoint has a 

143 transparent relationship to the local radar observations, whereas in a global solver, observations 

144 across the radar domain have an impact on the solution at an individual grid point. Without such 

145 interference, possible errors in the radar measurements or retrieved winds are more easily 

146 identifiable and traceable. With a natural coordinate system, the coplane method is particularly 

147 useful for understanding the advantages and disadvantages of the HIWRAP scanning technique. 

148 This understanding is necessary for interpreting any Doppler analysis method used on the 

149 HIWRAP radar geometry. In this study, we apply the coplane technique to simulated radar data 

150 and actual radar data to demonstrate its effectiveness. We also compare this technique to the 

151 global optimization solutions and investigate their differences. 

152 Sections 2-3 describe the coplane method and its application to the HIWRAP geometry 

153 and observations. Section 4 examines the boundary conditions necessary for the coplane method. 

154 Section 5 analyzes the coplane retrieval of simulated radar data and Section 6 analyzes the 

155 coplane retrieval of real HIWRAP data. Finally, Section 7 presents the conclusions of this study. 

156 

157 2. Coplane method and HIWRAP geometry 

158 a. Description of HIWRAP 

159 HIWRAP is a dual-beam, dual-frequency (Ka and Ku band) radar system designed to fly 

160 on the high-altitude NASA Global Hawk unmanned aircraft system. Rather than scanning around 

161 a horizontal axis like the tail radars on the P3 aircraft, the antenna beams of HIWRAP point 
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162 downward and scan around a vertical axis to obtain multiple angled looks of the tropospheric 

163 winds. Figure 1 illustrates this scanning geometry. The two beams point at nominal tilt angles (x) 

164 of 30° and 40° away from nadir, while the antenna rotates at a typical rate of 100° per second, 

165 such that one complete revolution takes about 3.5 s. The radar beams, each with a range 

166 resolution of 150 m, sweep out spiral paths on the ground as the aircraft flies with an ideal level 

167 position along a straight flight track. For a typical aircraft speed of 160 m s' 1 and altitude of 18.5 

168 km, the along-track sampling and swath width are 560 m and ~30 km respectively. The outer 

169 beam operates simultaneously at 13.5 and 33.7 GHz and the inner beam operates simultaneously 

170 at 13.9 and 35.6 GHz. HIWRAP employs dual pulse repetition frequency sampling that can yield 

171 an extended unambiguous velocity of ~1 10 m s' 1 . A more detailed description of HIWRAP can 

172 be found in Li et al. (2014). 

173 

174 b. Description of coplane method 

175 The coplane dual-Doppler technique was developed to retrieve the three-dimensional 

176 winds with two or more ground radars (Armijo 1969; Miller and Strauch 1974) and later applied 

177 to airborne tail radars employing the FAST scanning technique (Chong and Testud 1996). This 

178 technique is implemented in a cylindrical coordinate system whose central axis is the line 

179 between the location points where the radar (or radars) provides two different looks of a single 

180 point in the domain. For aircraft observations, the ideal situation for the coplane method would 

181 have a straight flight track and constant flight altitude across the analysis domain. The two looks 

182 of the wind field, obtained with fore and aft pointing beams, are considered independent and, for 

183 the purpose of this study, instantaneous. For the typical Global Hawk speed and altitude, the 

184 largest time gap between the observations is 200 s. These two measurements can then be readily 
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converted into two orthogonal wind components. Recovery of the third wind component at every 
point in the domain requires well-posed data (i.e., data exists at every point). To describe the 
application of the coplane method to the downward pointing conically scanning HIWRAP 
geometry, we follow the discussions from Tian et al. (2015) and Guimond et al. (2014). 

Figure 2 illustrates the cylindrical coordinate system defined by p, a, and Y. The flight 
track serves as the main axis Y where the origin is some arbitrary point along Y. The variable p is 
the radial distance from the central axis, and a is the coplane angle beginning at 0° for the nadir 
plane and increasing to the right of the flight track. For every rotation angle 9 (0° points in +Y 
direction), range r, and current track position Y 1 , observations are first mapped onto a track- 
following Cartesian grid by 

( X t \ / cosD(a) — sinD sinr(b) — sinD(c) \ 

T t ) = r( Y 1 /r + sinD(a) + cosD sinr(h) + cosD(c) J (1) 

Z t / V sin r (sin P cos 9 — cos P sin R sin 0) — cos P cos R cos r / 


where 


/ cl\ / cos R sin 9 sin t — sin R cos t \ 

( b J = I cos P cos 6 + sin P sin R sin 8 I (2) 

w' V sin P cos R cos r / 


and D, P, R and r are the drift, pitch, roll and tilt angles, respectively. Equations (1) and (2) are 
similar to Guimond et al. (2014) and are derived for the current scanning geometry following 


Lee et al. (1994). The cylindrical coordinates of the observations are then calculated by 



\ n / 


( 3 ) 


As the plane flies along the track, a single beam at a given tilt angle r obtains Doppler 
velocities in an a plane when it is located at Y 1 (fore) and Y 2 (aft). These velocities are 
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interpolated to the cylindrical coordinate grid so that each grid point P contains consolidated fore 
and aft radial velocities (17 rl , V r2 respectively) as seen in Fig. 2. Orthogonal velocities in the a 
plane are then calculated by 




-r 1 (Y-Y 2 )V rl +r 2 (Y-Y 1 )V r2 

p(X 2 - n) 

r l^rl — r 2^r2 


( 4 ) 


Uy 


y 2 -y 1 

where r x = p 2 + (Y — Y^ 2 and r 2 = -J p 2 + (Y — Y 2 ) 2 . From these standard dual-Doppler 
calculations, we obtain two velocity components (U p , U Y ) in each a plane of the cylindrical grid. 

The separation angle /?, defined as /? = /? x + /? 2 as seen in Fig. 2, is the angle between the 
fore and the aft beams. The angles /? x and /? 2 are calculated by 

p n = sin -1 (|y- Y n \/r n ) (5) 

where /? n represents either /? x or /? 2 . Combining Eqs. ( l)-(3) and assuming all attitude angles are 
equal to 0, Eq. (5) can be rewritten as 

P n — sin -1 (cos r e cos (sin -1 (tan r e tana))) (6) 

where /? n is now a function of the coplane angle a. Eq. (6) uses the elevation angle r e which is 
defined as r e — t — 90°. The separation angle directly corresponds to the accuracies of the two 
retrieved wind components, U p and U Y . In applying the error estimates of Doviak et al. (1976) 
and trigonometric substitutions to Eq. (4), the variances of the two wind components are 
specified by 


2 = g rl 2 + <b-2 2 2 = ^rl 2 + gr2 2 

° p Acos 2 ^ ° Y dsin 2 /?! 

where o rl and o r2 are the errors of V rl and V r2 . The errors o rl and o r2 are equal to each other 
given that V rl and V r2 are independent measurements. Tian et al. (2015) determined that the 
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217 standard error of HIWRAP Doppler estimates for the Ka band is a r = 0.46 ms" 1 . Eq. (7) assumes 

218 that all errors are Gaussian distributed. Other sources of error can contribute to o r such as 

219 velocity unfolding error and error due to aircraft motion. For the analysis in Section 6, we 

220 verified that the Doppler velocities were unfolded properly. The Doppler velocities were also 

221 corrected for aircraft motion using attitude infonnation (i.e. roll, drift, and pitch). 

222 Figure 3 shows the separation angle and the corresponding wind variances as a function 

223 of the coplane angle for the two tilt angles of the HIWRAP geometry. It is shown that fi reaches 

224 its peak at nadir and then decreases as a increases in magnitude. a Y 2 is lowest at nadir and 

225 remains below 0.6 m 2 s" 2 throughout most of the domain. Towards the domain edges, the fore and 

226 aft beams become closely parallel (i.e. approaches 0°) and point less in the along-track 

227 direction. Consequently, the accuracy of the retrieved U Y quickly degrades at large a magnitudes. 

228 On the other hand, the U p component is accurately estimated {o 2 < 0.2 m 2 s' 2 ). It is most 

229 accurate near the domain edges and least accurate at nadir. Still, the magnitudes of o p 2 and its 

230 corresponding changes with a are lower than that of o Y 2 . Studies have shown that the two in- 

231 plane wind components can both be retrieved with reasonable accuracy when the separation 

232 angle is at least 30° (e.g. Klimowski and Marwitz 1992). In this scanning geometry, the outer 

233 beam retrieves the wind components with reasonable accuracy when |a| < 37.5° where 

234 er 2 < 1.56 m s -1 for both components. 

235 Figure 3 also shows that the outer beam retrieves U Y more accurately while the inner 

236 beam retrieves U p more accurately within its smaller domain. We incorporate observations from 

237 both beams by weighting these relative retrieval accuracies. For each gridpoint within the 

238 domain of the inner beam, the composite wind components are 
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( 8 ) 


239 

240 

241 

242 

243 

244 

245 
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247 

248 

249 
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Un = 


®pi Upo 3 ” Opo Upi 


P q .2 -|_ (j 2 

u pi ' u po 


Uy = 


Oyi^Uvr, + Gyr^Uyl 


J Yi u YO T u Vo U yj 
a Yi 2 + °Yo 2 


where the i and o subscripts denote observations from the inner and outer beams. 

The third component of the wind, U a , is retrieved by integrating the anelastic mass 
continuity equation along the a-axis away from the nadir plane. Figure 4 illustrates the two 
integration directions that span the radar domain. The anelastic mass continuity equation is given 
by 

d(pvU r ) . d(rjU a ) d(rjU Y ) 


+ 


+ r- 


0 


dp da dY 

where >] is the air density. The current calculations use the Jordan (1958) standard Tropical 
Atlantic air density profile. Using the square rule for integration on Eq. (9), U a is obtained by 


( 9 ) 


V U a 


a , 


= VUa 


Uai - a 0 )(/(«i) +/Oo)) 


d{ppU r ) d(rjU Y ) 
/(a) = — Yr- 


( 10 ) 


dp dY 

where the subscripts 0 and 1 denote the previous and current integration locations. As depicted in 
Fig. 4, U a must be initialized with boundary conditions at the nadir plane and at the surface. In 
order to retrieve U a at all points, the data must exist at all points in the domain. If radial 
velocities are missing at any point, U a cannot be calculated at points along the integration path 
beyond the missing point. Data may continue beyond the missing point allowing for calculation 
of U Y and U p . 

The nadir boundary condition is obtained by taking observations at small angles away 
from nadir on either side. In the track-following Cartesian grid, the Cartesian coordinate cross- 
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track («), along-track (v), and vertical (w) velocities are related to the cylindrical coordinate 
velocities by 

U p — u sin a — w cos a (11) 

U Y = v ( 12 ) 

U a — u cos a + w sin a (13) 

Suppose that U p components (U pl , U p2 ) are calculated at a small angle a on either side of nadir 
(a 1 — +a; a 2 — —a) at a constant radius. For the two U p components, we make the assumption 
that w is constant and u is linear across the span of the U p locations. It follows from Eq. (11) that 
it at nadir (it 0 ) is expressed by 


Un 


Up 1 - Up 2 


(14) 


2 sin a 

for each altitude corresponding to the radius of the U p observations. Since U a = u 0 at nadir, Eq. 
(14) gives the boundary condition for initializing U a along the nadir plane. To calculate the nadir 
boundary condition, we chose a value of a = 3.35° for the outer beam and a = 2.31° for the 
inner beam. The two it 0 values from each beam at each point are combined according to the U p 
calculation and error estimates from Eq. (8). This weighting was selected since the final values 
rely on U p calculations. 

At the surface, the impermeability condition (w = 0) is applied as a boundary condition. 
By setting w equal to 0, Eqs. (11) and (13) lead to 


U n 


z = 0 


— U 0 / tan a 


(15) 


With this relationship, U p can be used to initialize U a at the surface. The surface boundary 
condition works well in an idealized setting where accurate observations are available near the 
surface and the surface is flat. However, in actual aircraft observations over water, sea spray can 
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271 contaminate the Doppler measurements and the surface is not flat. We address these surface 

272 issues and assess the nadir boundary condition in Section 4. 

273 

274 3. Data and methods 

275 a. Radar simulator 

276 In order to assess the validity of the coplane analysis, we use model output and a radar 

277 simulator designed after Guimond et al. (2014) with no added noise or aircraft attitude. The radar 

278 simulator mimics the scanning technique of the HIWRAP radar and obtains radial velocities 14- 

279 from the modeled velocity fields as the radar moves along a straight level track. The model used 

280 is the nonhydro static fifth-generation Penn State University- National Center for Atmospheric 

281 Research Mesoscale Model (MM5). We take an MM5 simulation of Hurricane Rita (2005) at a 

282 single time frame near its peak intensity (maximum wind speed of 75 m s' 1 ). The model output 

283 has a horizontal resolution of 1 .67 km and 28 sigma levels in the vertical. Two simulated radar 

284 beams are positioned at 30° and 40° tilt angles, and rotate at a period of 3.5 seconds per 

285 revolution with an azimuthal resolution of 2° and a range resolution of 150 m. The radar has a 

286 nominal altitude of 18.5 km and the simulated aircraft has a ground speed of 160 m s' 1 . Shown in 

287 Fig. 5, the track has a length of 200 km and passes through the center of the storm. 

288 

289 b. Real data 

290 On 16 September 2013, the NASA Global Hawk AV-1 flew over Hurricane Ingrid as part of 

291 the NASA Hurricane and Severe Storm Sentinel (HS3) field campaign. The HIWRAP radar 

292 on board the Global Hawk observed the northern edge of Ingrid as the stonn tracked west 

293 across the Gulf of Mexico. The data used in this study were taken from 1836-1900Z. Figure 6 
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294 


shows the HIWRAP observed reflectivity (plan view and at nadir) along with the 

295 corresponding infrared satellite image. In Section 6, we apply the coplane analysis to the Ka- 

296 band outer beam observations. To remove noise, pixels with reflectivity less than 0 dBZ were 

297 not used in the analysis. The Doppler velocities were unfolded according to Dazhang et al. 

298 (1984). We applied corrections for beam pointing errors by aligning the expected range of 

299 the ocean surface with the range of the observed surface return. Fall speed corrections from 

300 Heymsfield et al. (2010) were also applied to the velocity data. In this correction algorithm, 

301 fall speeds were calculated as a function of the Ka band reflectivity and altitude. 

302 c. Grid and interpolation specifications 

303 The coplane method requires an initial interpolation of radial velocity data to a 

304 cylindrical grid. The cylindrical grid used in this study has a radial resolution of 0.5 km, along- 

305 track resolution of 2 km, and azimuthal resolution of 2.5°. The observations (both simulated and 

306 real) are interpolated to this grid using a Bames weighting scheme (Bames 1973; Koch et al. 

307 1983) given by 

w m - exp ) (16) 

308 where r m is the distance of the m th observation from the analysis grid point, y is a chosen shape 

309 parameter, and 6 is the influence radii expressed by 



310 where r p ,r Y , and r a are the radii of influence in the three coordinate directions. For this 

311 interpolation, the radial, horizontal, and azimuthal radii of influence are 0.5 km, 2 km, and 1.25° 

312 respectively. The 1 .25° azimuthal radius of influence has an equivalent distance of 2 km and the 

313 shape parameter y , which determines the width of the weighting function, is chosen as 0.75. 
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314 Following the coplane calculations, the data are converted into Cartesian coordinate velocities 

315 via Eqs. ( 1 1)-( 13) and are finally interpolated to Cartesian coordinates. The Cartesian grid has a 

316 horizontal resolution of 2 km and a vertical resolution of 1 km. An additional level is added at 

317 0.5 km altitude for better resolving of the low-level winds. This interpolation uses the same 

318 Barnes filter but with a radius of influence of 2 km in the horizontal dimensions and 0.25 km in 

319 the vertical in the vertical dimension. By determining the response function of the Bames filter 

320 (Koch et al. 1983), the minimum resolvable horizontal wavelength is calculated to be 4 km, 

321 which is also twice the horizontal grid spacing. 

322 

323 4. Boundary conditions analysis 

324 a. Nadir boundary condition analysis 

325 At nadir, the cross-track wind component is unobserved by the HI WRAP radar and must 

326 be estimated by utilizing other available measurements. It is important to obtain a good estimate 

327 of the cross-track wind at nadir as this will serve as the boundary condition that initializes the U a 

328 wind component for integration throughout most of the domain. As shown in Eq. (14), we 

329 estimate the cross-track wind by using wind measurements taken at a small angle a away from 

330 nadir. Choosing a value for a requires a balance of certain trade-offs. For smaller a values, the 

331 distance between observations is smaller and thus the assumptions of constant vertical velocity 

332 and linear cross-track velocity are well suited. However, at angles that are closer to zero, the 

333 wind measurements are more susceptible to errors in the cross-track velocity. For larger a 

334 values, the cross-track wind is better sampled and this reduces the susceptibility to measurement 

335 errors; however, the distance between observations is greater making the necessary assumptions 

336 less suitable. 


15 



337 We use the simulated radar data to choose a value for a. Boundary condition estimates 

338 are calculated with varying a values, which are then compared to the model “truth” cross-track 

339 velocities. This calculation requires interpolation of radial velocities to the different a planes. A 

340 Barnes filter is used for the interpolation with the influence radii specified in Section 3c. 

341 Assuming all attitude angles are equal to zero, Eqs. (l)-(3) yield a as a function of the rotation 

342 angle 9 and the elevation angle r e : 

„ /sin 9 cos r„\ 

a = tan -1 ( : - (18) 

V sin r e ) 

343 The a values are tested by varying the deviation of 9 from the nadir plane. For example, the 

344 rotation angles 2° and 178° lie in the plane a = 1.677° for the outer beam. Correspondingly, the 

345 rotation angles 358° and 182° lie in the plane a — —1.677°. Figure 7a shows the coplane angles 

346 for the varying rotation angle, while Figure 7b presents the root-mean-square (RMS) errors for 

347 the different estimates. The errors are all relatively small compared to the wind speeds of the 

348 simulated hurricane. For a wind speed of 30 m s' 1 , the largest error in Fig. 7b constitutes 5% of 

349 this wind speed. The a values corresponding to 9 = 4° produced the smallest error, so we 

350 chose these values for the boundary condition retrieval. As seen in Fig. 7b, the angle 9 — 4° 

351 corresponds to a cross-track distance between observations of 2.2 km at the surface. 

352 Figure 8 displays the estimated U a at nadir along with the errors relative to the model 

353 truth. The retrieved wind field captures the overall structure of the hurricane. Errors larger than 2 

354 ms' 1 occur near the eyewall region (7= 85 and 115 km), the surface, and the domain edge at 7 = 

355 18 km. These positive errors at the domain edge reach 8 m s' 1 . The largest negative errors occur 

356 in the midlevels near 7= 155 km reaching values of -9 ms' 1 . These errors stem from local 

357 violations of the assumptions made in the calculation of Eq. (14). Specifically, violations of the 
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358 constant vertical velocity assumption are the primary source of error in Fig. 8, where vertical 

359 velocity deviations of 0.5 m s' 1 produced U a errors of ~3 ms' 1 . 

360 

361 b. Surface boundary condition analysis 

362 The values of U a must be initialized at the lower boundary of the analysis domain. Given 

363 the curved paths of integration, this initialization affects the lower portion of the domain that 

364 increases in depth away from nadir (as indicated in Fig. 4). As described in Section 2, the lower 

365 boundary initialization can be done most simply by invoking the impermeability condition and 

366 setting w = 0 at the surface (Eq. 15). This condition requires reliable observations near a flat 

367 surface, which is an ideal situation that models provide. The initialization locations on the 

368 surface are not necessarily points on the cylindrical grid, but U a can still be effectively initialized 

369 for every integration path that intersects the surface. 

370 With actual observations, setting the surface boundary condition cannot be done so 

371 simply, particularly over the ocean surface as in the case of tropical cyclone research. The ocean 

372 surface may not be flat and sea spray can contaminate echoes near the surface. Previous dual- 

373 Doppler methods approach the surface initialization of the integrated wind component (usually 

374 w) differently. In the Cartesian Editing and Display of Radar Data under Interactive Control 

375 software (CEDRIC; Mohr et al. 1986), the vertical velocity can be initialized at the lowest level 

376 of usable data by setting w equal to a specified multiple of the locally measured horizontal 

377 divergence. Chong and Testud (1996) use a variational method in which w at all surface 

378 locations vary near w = 0 such that the resulting U a field contains a minimum amount of 

379 continuity irregularities. 
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For this study, we use a simple approach that initializes U a at the lowest cylindrical grid 
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points closest to a selected low-level altitude. In the simulated data, we choose 0.5 km as the 
lowest level of available data, which is approximately the lowest level of usable data from the 
HIWRAP observations. For this cylindrical coordinate system, the lower boundary grid points 
are not all at the same altitude. By combining Eqs. (11) and (13), U a at each lower boundary 
point is given by 


U p ( cos* a 

J„ — h w 1- sin a 

tana \ sin a 


Using this equation, U a is initialized with the local U p . Additionally, we estimate w from the 
vertical velocity calculated at nadir (where w = —U p ) for that corresponding altitude and Y 
location. 


The calculated lower bound U a values were compared to the model truth and resulted in 
an RMS error of 2.05 m s' 1 . Moreover, the lower bound U a values stemming from the original 
impermeability condition resulted in an RMS error of 1.69 m s' 1 . As expected, the accuracy of 
the nadir-vv approach is slightly lower than that of the impermeability approach; however, the 
difference in the errors (0.36 m s' 1 ) is small relative to the near-surface hurricane wind speeds 
(which have magnitudes greater than 30 m s' 1 outside the eye. From this analysis, the nadir-vv 
approach represented by Eq. (19) is deemed suitable for U a initialization for near-surface grid 
points. 


5. Retrieval error analysis 

In this section, we use the simulated radar data to examine the wind field retrieved from 
the coplane analysis. As described in the previous section, we do not use radar radial velocities 
below 0.5 km altitude in this retrieval. Figures 9a-c present the RMS errors calculated along the 
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402 flight track for the retrieved cross-track (u), along-track (v), and vertical velocity (w) 

403 components. These figures show the total errors and the error patterns of each wind component 

404 for the HIWRAP scanning geometry. The total relative-RMS (RRMS) error in each figure is a 

405 nonnalization of the errors relative to the velocity magnitudes. 

406 The cross-track component (u) contains an average error of 1.9 m s' 1 which, as indicated 

407 by the relative-RMS value of 4.4%, is a low error compared to the u magnitudes. Calculation of 

408 u depends on both the U a and U p components, but the u errors largely stem from errors in U a as 

409 this component is larger and more aligned with u throughout the domain. The u errors form a 

410 curved pattern as they follow the integration path upon which U a was calculated. The largest 

411 errors occur near the surface and in a midlevel belt positioned between 4-6 km altitude at nadir. 

412 The vertical velocity (w) contains an average error of 0.9 m s' 1 which is significant relative to the 

413 vertical velocity magnitudes (RRMS = 60.4%). Despite this significant average error, the error 

414 distribution in Fig. 9c shows that the vertical velocities near nadir have the smallest errors and 

415 therefore are the most useful. The errors increase as the a angle magnitude increases toward the 

416 edges of the domain, with particularly large magnitudes at locations that coincide with the u error 

417 belt in Fig. 9a. At these larger a angles, U a makes an increasing contribution to determining w. 

418 Asa result, U a errors that are small relative to the horizontal winds can lead to significant w 

419 errors near the domain edges. 

420 We have explained that errors in the u and w fields are mostly due to U a errors. These 

421 errors in the U a component accumulate during the integration of the wind field for two reasons. 

422 First, the divergence of the wind field in the a planes is not well sampled, particularly near the 

423 domain edges where U Y calculations become less accurate (Fig 3b). Second, U a is incorrectly 

424 initialized for the two boundary conditions. We briefly test which reason is most responsible for 
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425 the U a errors by substituting the lower bound and nadir U a estimates with the model truth. 

426 Figure 10 shows the RMS error patterns. Having the best initialization possible, the wind field 

427 errors are significantly reduced to 1.1 m s' 1 for u and 0.5 m s' 1 for w. The it error no longer 

428 contains the belted pattern and the corresponding w errors along the domain edges are removed. 

429 This analysis suggests that the errors in the u and w wind fields are mostly a result of errors in 

430 the boundary conditions. The remaining errors are less pronounced in Fig. 10 and can be 

431 attributed to divergence sampling and interpolation error. 

432 The along-track component (v) is the only Cartesian coordinate component that is not 

433 calculated with the U a component. The v errors are very small throughout most of the domain. 

434 The largest errors, reaching up to 5 m s' 1 , occur at 0.5 km altitude. These errors at the lowest 

435 level are largely a result of the interpolation from cylindrical to Cartesian coordinates. The lower 

436 bound points on the cylindrical grid are at different altitudes and are all higher than the 0.5 km 

437 level. Since the data below 0.5 km were not used, these lower bound points are the only source 

438 of information for interpolation onto the 0.5 km level, which contributes to errors found at this 

439 lowest level for all Cartesian wind components. In this particular dataset, the v component 

440 (which is largely the radial wind of the hurricane) changes very rapidly at these boundary layer 

441 altitudes, resulting in the significant errors found in the lowest levels. When the 0.5-km-level 

442 data are excluded from the error analysis, the RMS error drops from 1.7 m s' 1 to 1.0 m s' 1 . 

443 The coplane analysis is now compared to the global optimization analysis described by 

444 Guimond et al. (2014). In their variational method, radar velocities are first interpolated to a 

445 Cartesian grid, and then a modeled wind field is retrieved using the radar velocities, mass 

446 continuity, and boundary conditions as constraints. We applied the variational method to the 
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447 current simulated data using analysis parameters that produced the smallest errors. Specific 

448 parameters that were chosen are listed in Table 1 . 

449 Figures 9d-f show the RMS errors for the three wind components. The u and w 

450 components both have larger overall errors than the coplane analysis. These components, which 

451 again are connected to each other through the unobserved wind component U a , have error 

452 patterns that are slightly different from the coplane analysis errors. The best retrievals occur at 

453 nadir, and errors increase at all altitudes when approaching the domain edges. These error 

454 patterns are consistent with those from Guimond et al. (2014). There is no belt of errors as in Fig. 

455 9a, but increased errors do exist at nadir at the same altitude range of 4-6 km. 

456 Meanwhile, the v field has a smaller overall error than the coplane analysis. The largest 

457 difference in the v error pattern is in the lowest levels. The variational method operates fully in a 

458 Cartesian coordinate system, which means that the lowest level of available data coincides with 

459 the lowest level of grid points at 0.5 km altitude. As a result, the rapid vertical changes in v are 

460 well captured at these levels and there is no interpolation error from changing coordinate 

461 systems. 

462 

463 6. Coplane retrieval with real radar data 

464 In the previous section, the coplane retrieval method was successfully applied to 

465 simulated radar data. We now apply the coplane method to real HIWRAP data shown in Fig. 6 

466 and we compare the retrieved wind field to a solution from the variational method (Guimond et 

467 al. 2014). 

468 Figures 1 la-c presents the coplane analysis cross-track (u), along-track (v), and vertical 

469 (w) components of the wind field along nadir of the observation domain. The u field contains 
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470 mostly positive velocities with values > 10 m s' 1 in the layer below 3 km and the layer above 6 

471 km. In between these layers a midlevel minimum of u occurs. The v field also contains mostly 

472 positive velocities that decrease towards the upper levels of the domain (>7 km altitude). When 

473 considering the domain location (Fig. 6a), the u and v fields show consistency with the counter- 

474 clockwise cyclonic circulation of the storm. The cross-track component (Fig. 1 la) largely 

475 switches from positive to negative values at around F,=230 km which corresponds to the point 

476 along the track that is closest to the storm center. Concurrently, the along-track component (Fig. 

477 lib) increases as the track approaches the same closest point. 

478 The reflectivity field shown in Fig. 6c contains a clear bright band signature (at -4.5 km 

479 altitude) and fall streaks, which indicate that the dominant precipitation regime for these 

480 observations is stratiform (Houze 1997). In stratiform precipitation, falling ice crystals melt in a 

481 layer beneath the 0° isotherm and form a broad region of light to moderate precipitation. The w 

482 field (Fig lie) shows consistent features with stratifonn precipitation, including small 

483 magnitudes (< 2 m s' 1 ) throughout most of the domain. Updrafts are dominant above the bright 

484 band. Below this level, downdrafts are prominent, but a clear exception of positive w values 

485 occurs toward the beginning of the domain and at 4 km altitude. These exceptions, which are 

486 inconsistent with typical stratiform kinematics, are likely a result of errors in the fall speed 

487 correction and/or attenuation of the Ka beam. If these errors were consistent across the radial 

488 velocities used in Eqs. (4) and (14) to calculate U Y and u 0 , then the errors would not have an 

489 impact on the fields in Figs. 11a and lib. 

490 Figures 12a-c present the wind components along a cross section at Y t = 160 km. These 

491 cross sections show that the overall patterns seen at nadir extend to the edges of the domain. The 

492 midlevel minimum of u grows larger to the right of the flight track (Fig. 12a). Additionally, the 
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493 downdraft layer (Fig. 12c) also increases in depth to the right of the flight track. The slanted 

494 stretches of downdrafts < -3 m s' 1 near the domain edges do not appear consistent with expected 

495 vertical velocity patterns of stratiform precipitation. Rather, these patterns are reminiscent of the 

496 curved error patterns in Figs. 9a and 9c. Given this resemblance, we infer that these w swaths 

497 (and their corresponding u values) contain errors for the same reasons as in the radar simulator 

498 analysis - that is, incorrect initializations of U a at the boundary condition. By following the 

499 curved a paths from these features towards the domain center, one finds that the boundary 

500 condition errors occur at nadir between 1-2.5 km altitude in this cross section. 

501 The variational method retrieved a qualitatively similar wind solution as that of the 

502 coplane method. Figures 1 ld-f show that the u and v fields in the nadir plane have the same 

503 overall structures as in Fig. 1 la-c. Figure 1 If mostly has the vertical velocities expected of 

504 stratiform precipitation, but this solution takes the same fall speed corrections as in the coplane 

505 method and produces noticeably different vertical velocities at the bright band altitude (-4.5 

506 km). The variational u field (Fig. 1 id) has noticeably smoother contour patterns than the coplane 

507 u field (Fig. 11a), which suggests that the variational method may be filtering out some small- 

508 scale features in the data. The v and w fields from both methods do not have a noticeable 

509 discrepancy in their contour smoothness. Upon closer inspection, the it field at nadir is impacted 

510 most by the smoothing parameter in the variational retrieval. When this smoothing parameter is 

511 turned off, the resulting u field appears very similar to the coplane u field. 

512 One advantage of the coplane method is that the minimum resolvable wavelength of the 

513 data field is readily determined by calculating the response function of the Bames filter. On the 

514 other hand, determining the minimum resolvable wavelength of the variational method solution 

515 is not as straight forward. While a Bames filter is also used, the weighting parameter is a 
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516 constraint on the optimization and not a direct calculation. Thus the Barnes filter response 

517 function cannot exactly detennine the minimum resolvable wavelength. In addition, the 

518 smoothing parameter certainly increases the minimum resolvable wavelength, but again, this 

519 smoothing is a constraint and not a direct calculation. For both the Barnes filter and the 

520 Laplacian smoother, the minimum resolvable wavelength must be determined empirically. 

521 Figures 12d-f display the same cross section as in Figs. 12a-c but for the variational 

522 solution. As in the coplane analysis, the midlevel u minimum and downdraft layer increase in 

523 depth to the right of the flight track. The w field in Fig. 12f does not contain the unrealistic 

524 downdraft patterns seen in Fig. 12c as there is no explicit integration along a curved path. 

525 However, the w field does contain downdrafts < -3 m s' 1 near the domain edges that appear 

526 unrealistic. Along the left domain edge, these increased downdrafts occur in the same location as 

527 in the coplane analysis (Fig. 12c). Along the right domain edge, these increased downdrafts are 

528 prominent in the lower altitudes and appear to trail off into the higher altitudes. This pattern of 

529 vertically oriented anomalies along the domain edge is reminiscent of the error pattern in Fig. 9f, 

530 which suggest that these features contain likely errors. The source of these errors cannot be 

531 traced to specific observations, but rather, the errors must be attributed to the general decreased 

532 accuracy of the global solver along the domain edges. 

533 Both the coplane and variational methods produced adequate wind fields that generally 

534 agreed well with each other. Both fields also contained inevitable localized errors. With a priori 

535 knowledge of the error patterns expected from each method, the questionable features that appear 

536 in the solutions can be easily identified as retrieval errors. Identifying and understanding these 

537 errors is essential for reliable scientific interpretations of solutions from either analysis method. 
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538 We make a final comparison of retrieval techniques to the VAD technique from Lin et at. 

539 (2014). The VAD technique obtains the mean horizontal wind within the nadir plane by fitting 

540 the measured radial winds in a scan circle to a sinusoidal curve. Figure 13 shows the retrieved u 

541 and v components of the wind for the same leg of data from Hurricane Ingrid. The VAD 

542 technique captures the same overall wind pattern that was retrieved by the other retrieval 

543 techniques (Fig. 1 1). The most noticeable difference in Fig. 13 is the increased vertical 

544 resolution. Since the VAD technique does not retrieve the full three-dimensional wind field, it is 

545 computationally less expensive than both the coplane and variational methods; moreover, this 

546 allows the VAD technique to preserve the high vertical resolution of the HIWRAP beam. 

547 In calculating the mean horizontal wind, the wind field is assumed to have linear horizontal 

548 velocity and constant hydrometeor vertical speed across the total scan circle. These assumptions 

549 tend to hold well in stratiform precipitation regions like that in the current dataset since these 

550 mesoscale regions contain weak vertical velocities and winds that vary slowly over horizontal 

551 distances. In order to capture convective-scale features, one of the three-dimensional wind 

552 retrieval methods must be used. 

553 

554 7. Conclusions 

555 In this paper, the coplane method for dual-Doppler wind retrieval (Armijo 1969; Miller 

556 and Strauch 1974) is adapted to the downward pointing conically scanning technique of the 

557 NASA HIWRAP airborne radar. The coplane method takes the radar observations and solves for 

558 the three-dimensional winds using a simple interpolation and integration approach. This 

559 approach locally solves for the wind field which is in contrast to the global optimization 

560 (variational) method described by Guimond et al. (2014). In order to retrieve the unobserved 
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561 wind component ( U a ) at all points, observations must exist at all points in the domain. The main 

562 advantage of the coplane method is the transparency of its calculations. The interpolation and 

563 solving processes are discretely and separately calculated, which allows for exact calculation of 

564 wavelength resolution and tracing of source data from the solution. 

565 Simulated radar observations of a model hurricane were used to test the coplane method 

566 and compare to the variational method. The coplane method retrieved the wind field with small 

567 errors relative to the wind speed magnitudes. Compared to the variational method, the coplane 

568 method had lower errors in the cross-track component ( u ) and vertical component (w) fields, 

569 while the variational method had lower errors in the along-track component (v) field. For the 

570 coplane method, the accuracy of u relied on the accuracy of the U a boundary initializations. 

571 Where U a was initialized sufficiently well, u remained accurate across the span of the domain. 

572 Where U a was not well initialized, errors in u propagate along the curved integration path, 

573 creating an easily recognizable error signature. The w component, which is also derived from U a , 

574 produced errors at the domain edges along curved integration paths with insufficiently initialized 

575 U a . The error patterns for the variational field were different, showing errors in u and w that 

576 grew toward the edges of the domain at all altitudes. 

577 The coplane and variational methods were applied to HI WRAP observations collected 

578 during the NASA HS3 campaign. Both techniques produced errors in the retrieval that appeared 

579 in patterns similar to the errors in the simulated radar retrieval. Prior knowledge of the error 

580 patterns expected from each method allowed for this recognition of retrieval errors in the HS3 

581 retrieval. As a local solver, the errors in the coplane analysis are easily traced to the certain 

582 observations and/or U a initializations. Unlike the transparency of a local solver, the errors 

583 arrived with the global solver cannot be explicitly traced to certain observations or calculations 
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584 since the solution at a particular location depends on the solution everywhere. Additionally, since 

585 the coplane method employs the Bames filter, the corresponding response function provides the 

586 exact minimum resolvable wavelength of the final solution. The minimum resolvable wavelength 

587 for the variational method cannot be exactly calculated since the interpolation filter and 

588 Laplacian smoother are constraints on the optimization rather than exact calculations. 

589 The coplane technique’s ability to transparently trace the exact calculations from the raw 

590 observations to the final solution is highly beneficial when making scientific interpretations. This 

591 ability is necessary to more adequately analyze small-scale features in tropical cyclones, such as 

592 rotating deep convection (Hendricks et al. 2004; Montgomery et al. 2006; Sanger et al. 2014). A 

593 key skill for radar analyses is being able to separate true meteorological signals from non- 

594 meteorological signals such as noise or data contamination that has bypassed the data quality 

595 control process. Once a solution is obtained, quirky regions in the final solution can be 

596 ambiguous as to whether they are true representations of small-scale features. After easily 

597 pinpointing the raw observations that were used to create the solution, the user can better assess 

598 the reality of the observations taken, and corrections can subsequently be implemented or not 

599 implemented. In the variational method, these anomalous measurements would be smoothed and 

600 impact the entire retrieval, which would either dilute the small-scale signal or incorporate 

601 erroneous data in the solution. The coplane method can prevent incorrect scientific 

602 interpretations of inherently wrong solutions or strengthen confidence in the conclusions based 

603 on the observations. Given the wind component error analyses for both retrieval methods, the 

604 option to trace solution calculations with the coplane analysis is provided to the user with 

605 minimal cost to the accuracy of the overall solution. 
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606 Future work will use the coplane analysis for scientific research of observations from the 

607 HS3 field campaign and other planned campaigns using the HIWRAP radar. This technique can 

608 also be used to analyze tail Doppler radar data from tropical cyclones documented in peer- 

609 reviewed work (e.g. Marks and Houze 1984; Reasor et al. 2009; Houze et al. 2009; Bell and 

610 Montgomery 2010) as well as future field campaigns. In locations directly beneath the aircraft, 

611 the geometry of the tail Doppler radar observations is compatible with the coplane technique. 

612 These additional analyses would be especially useful here as these locations are particularly 

613 troublesome for capturing small-scale features with the global optimization technique. 

614 
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Tables 


720 

721 Table 1 . Parameters used for the variational dual-Doppler retrieval. Each parameter is explained 

722 in detail in Guimond et al. (2014). 

723 


Shape parameter y 

0.75 

Along-track sampling s 

560 m 

Smoothing factor /? 

4 

Weighting factor cxm 

2x2000 2 

Weighting factor as 

For simulated 
data: 0 

For real data: 
0.4x2000 4 
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739 Figure Captions 

740 

741 Figure 1. Schematic showing the scanning technique of the HIWRAP radar in a track-following 

742 Cartesian coordinate frame (X t , Y t , Z,). Two beams each with Ku and Ka band point downward at 

743 two angles and scan conically around a vertical axis. For a stationary radar and plane at its 

744 typical altitude of 18.5 km, the outermost beam scans a circle at the surface (Z=0) with an 

745 approximate diameter of 30 km. 

746 

747 Figure 2. Cylindrical coordinate system used for the coplane dual-Doppler retrieval from 

748 HIWRAP observations. P is the position of the target in the cylindrical coordinate system 

749 defined by the radius (p), coplane angle (a), and flight track distance (7). The red dashed lines 

750 and red arrows represent line segments and vectors, respectively, along the radar beam. The 

751 black dashed line and black arrows within the coplane represent a line segment and vectors in the 

752 cylindrical coordinate system. The symbols r\ and r 2 are the ranges of the target from the radar 

753 positions at Y\ and Yi. f>\ and /fi are the positive angles between the ranges and the radius p. V x \ 

754 and Vti are the radial Doppler velocities at point P, and U p and Uy are the corresponding 

755 cylindrical coordinate velocity components in the coplane. The axes X t , Y, and Z t define the 

756 corresponding track- following Cartesian coordinate frame where T t is coincident with Y. Figure 

757 adapted from Tian et al . (2015). 

758 

759 Figure 3. a) The separation angle /? as a function of the coplane angle a for the inner beam 

760 (dashed black) and the outer beam (solid black), b) Variance (cr 2 ) of the Uy (red) and U p (blue) 

761 wind components due to Doppler estimate error as a function of the coplane angle a. Variances 


34 



762 are calculated for the Ka band inner (red dashed and blue dashed) and outer (red solid and blue 

763 solid) beams. 

764 

765 Figure 4. Diagram of the cylindrical grid ( p , a ) along a given 7 location. Red lines indicate the 

766 boundary condition locations for initialization of the U a wind component. Blue arrows indicate 

767 the integration directions for retrieving the U a wind component. The shaded region shows the 

768 area initialized by the surface boundary condition. 

769 

770 Figure 5. Plan view of model reflectivity at 2 km altitude. The data are HIWRAP radar simulator 

771 observations of a model simulation of Hurricane Rita at peak strength. 

772 

773 Figure 6. a) Infrared satellite imagery of Hurricane Ingrid (2014) with the flight track from 15 

774 Sep 2014 of the NASA Global Hawk AV 1 overlain in yellow. A black ‘X’ marks the best track 

775 storm center b) The 1 .5 -km Ka band outer beam reflectivity observed by HIWRAP along the 

776 flight track from a). Observations were taken from 1836-1900Z. c) Nadir view of Ka band outer 

777 beam reflectivity. The graph origin is the eastern end of the plan view reflectivity from b). 

778 

779 Figure 7. a) The coplane angle a as a function of rotation angle calculated from Eq. (18). 

780 Results from the inner beam (red) and outer beam (black) are shown, b) RMS errors of U a in the 

781 nadir boundary condition estimate for varying rotation angles (dots) and cross-track distance 

782 between outer beam locations for varying rotation angles (dashed line). 

783 
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784 Figure 8. U a wind component at nadir retrieved from the radar simulator data using the coplane 

785 analysis. Deviations from the model truth are shown in black contours at intervals of 2 m s' 1 . 

786 Dashed lines are negative values beginning at -2 m s' 1 and solid lines are positive values 

787 beginning at 2 m s' 1 . See text for details. 

788 

789 Figure 9. Root-mean-square (RMS) errors of the two retrieval methods. Errors from the coplane 

790 analysis are shown for the a) cross-track (u), b) along-track (v), and c) vertical (vv) wind 

791 components. Errors from the variational analysis are shown in d) - f) for the same wind 

792 components. Error fields are averaged along the flight track. Overall RMS and relative-RMS 

793 errors are also given (in m s' 1 and %, respectively). 

794 

795 Figure 10. RMS errors of the coplane analysis for the a) cross-track ( u ) and b) vertical (vv) wind 

796 components averaged along the flight track. This analysis uses nadir and lower-bound boundary 

797 conditions given by the model truth field. Overall RMS and relative-RMS errors are also given 

798 (in m s' 1 and %, respectively). 

799 

800 Figure 11. Nadir view of the a) u, b) v, and c) w wind components as derived by the coplane 

801 analysis of the HIWRAP observations seen in Fig. 6. The variational analysis wind components 

802 are shown in d) - f). 

803 

804 Figure 12. Cross-track view of the a) u, b) v, and c) vv wind components as derived by the 

805 coplane analysis of the HIWRAP observations. The variational analysis wind components are 

806 shown in d) - f). This cross section is taken at E = 160 km from Fig. 12. 
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807 


808 Figure 13. Nadir view of the a) u and c) v wind components as derived by the vertical azimuth 

809 display (VAD) analysis of the HIWRAP observations seen in Fig. 6. 
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Figures 


HIWRAP scanning geometry 



Figure 1 . Schematic showing the scanning technique of the HIWRAP radar in a track-following 
Cartesian coordinate frame {X,, Y t , Z,). Two beams each with Ku and Ka band point downward at 
two angles and scan conically around a vertical axis. For a stationary radar and plane at its 
typical altitude of 18.5 km, the outermost beam scans a circle at the surface (Z= 0) with an 
approximate diameter of 30 km. 
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Figure 2. Cylindrical coordinate system used for the coplane dual-Doppler retrieval from 
HIWRAP observations. P is the position of the target in the cylindrical coordinate system 
defined by the radius (p), coplane angle (a), and flight track distance (7). The red dashed lines 
and red arrows represent line segments and vectors, respectively, along the radar beam. The 
black dashed line and black arrows within the coplane represent a line segment and vectors in the 
cylindrical coordinate system. The symbols r\ and r-i are the ranges of the target from the radar 
positions at Y\ and Y 2 . f>\ and /1 2 are the positive angles between the ranges and the radius p. V r \ 
and V r 2 are the radial Doppler velocities at point P, and U p and Uy are the corresponding 
cylindrical coordinate velocity components in the coplane. The axes X t , Y, and Z t define the 
corresponding track-following Cartesian coordinate frame where T t is coincident with Y. Figure 
adapted from Tian et al. (2015). 
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Figure 3. a) The separation angle /? as a function of the coplane angle a for the inner beam 
(dashed black) and the outer beam (solid black), b) Variance (a 2 ) of the Uy (red) and U p (blue) 
wind components due to Doppler estimate error as a function of the coplane angle a. Variances 
are calculated for the Ka band inner (red dashed and blue dashed) and outer (red solid and blue 
solid) beams. 
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Cylindrical grid for 
coplane analysis 



Figure 4. Diagram of the cylindrical grid (p, a) along a given 7 location. Red lines indicate the 
boundary condition locations for initialization of the U a wind component. Blue arrows indicate 
the integration directions for retrieving the U a wind component. The shaded region shows the 
area initialized by the surface boundary condition. 
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Figure 5. Plan view of model reflectivity at 2 km altitude. The data are HIWRAP radar simulator 
observations of a model simulation of Hurricane Rita at peak strength. 
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Figure 6. a) Infrared satellite imagery of Hurricane Ingrid (2014) with the flight track from 15 
Sep 2014 of the NASA Global Hawk AV 1 overlain in yellow. A black ‘X’ marks the best track 
storm center b) The 1.5-km Ka band outer beam reflectivity observed by HIWRAP along the 
flight track from a). Observations were taken from 1836-1900Z. c) Nadir view of Ka band outer 
beam reflectivity. The graph origin is the eastern end of the plan view reflectivity from b). 
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Rotation angle 0 

Figure 7. a) The coplane angle a; as a function of rotation angle calculated from Eq. (18). 
Results from the inner beam (red) and outer beam (black) are shown, b) RMS errors of U a in the 
nadir boundary condition estimate for varying rotation angles (dots) and cross-track distance 
between outer beam locations for varying rotation angles (dashed line). 
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Figure 8. U a wind component at nadir retrieved from the radar simulator data using the coplane 
analysis. Deviations from the model truth are shown in black contours at intervals of 2 m s' 1 . 
Dashed lines are negative values beginning at -2 m s' 1 and solid lines are positive values 
beginning at 2 m s' 1 . See text for details. 
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Figure 9. Root-mean-square (RMS) errors of the two retrieval methods. Errors from the coplane 
analysis are shown for the a) cross-track (u), b) along-track (v), and c) vertical (vv) wind 
components. Errors from the variational analysis are shown in d) - f) for the same wind 
components. Error fields are averaged along the flight track. Overall RMS and relative-RMS 
errors are also given (in m s' 1 and %, respectively). 
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Figure 10. RMS errors of the coplane analysis for the a) cross-track ( u ) and b) vertical (vv) wind 
components averaged along the flight track. This analysis uses nadir and lower-bound boundary 
conditions given by the model truth field. Overall RMS and relative-RMS errors are also given 
(in m s ’ 1 and %, respectively). 
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935 

936 Figure 11. Nadir view of the a) u, b) v, and c) w wind components as derived by the coplane 

937 analysis of the HIWRAP observations seen in Fig. 6. The variational analysis wind components 

938 are shown in d) - f). 
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Figure 12. Cross-track view of the a) u, b) v, and c) w wind components as derived by the 
coplane analysis of the HIWRAP observations. The variational analysis wind components are 
shown in d) - f). This cross section is taken at Y t = 160 km from Fig. 12. 
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949 Figure 13. Nadir view of the a) u and c) v wind components as derived by the vertical azimuth 

950 display (VAD) analysis of the HIWRAP observations seen in Fig. 6. 
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